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■ Abstract 

A general fundamental relation connecting the correlation of Stokes and 
anti-Stokes modes to the quantum statistical behavior of vibration and pump 
modes in Raman-active materials is derived. We show that under certain 
conditions this relation can be used to determine the equilibrium number 
variance of phonons. Time and temperature ranges for which such conditions 
can be satisfied are studied and found to be available in todays' experimental 
standards. Furthermore, we examine the results in the presence of multi-mode 
\ pump as well as for the coupling of pump to the many vibration modes and 

O ' discuss their validity in these cases. 
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I. INTRODUCTION 



The concept of squeezed state has been established in the language of physics mainly 
by the developments in quantum optics. On the other hand, basic requirement of finding a 
system in a squeezed state is to have bosons as the constituents of the system interacting in 
a pairwise manner and that might be fulfilled not only in optical systems but in some other 
Bose-type systems as well. In actual fact, the introduction of squeezed states in opticsS 
was based on the previous consideration of superfluidityi in liquid He 4 (also see§). While 
squeezing of quantum fluctuations is the most well-known aspect of squeezed states, rich 
variety of effects might be expected due to their interesting statistical properties even at 
thermal equilibrium. Certain effects like anti-bunching have already been observed in the 
realm of quantum optics and this makes it an intriguing question how to find squeezed states 
and their effects in other places. In this context, few proposals have been suggested for the 
generation and detection of squeezed states of Bose-type excitations in solids@~i. Quite 
recently, squeezed phonons have been produced and detected^. 

It is very interesting that, unlike the case of light, the squeezed states of phonons may 
arise from different microscopic interactions in solids even at thermal equilibrium!. Devia- 
tions from typical equilibrium distribution of phonons, namely Bose-Einstein distribution, 
might arise from anharmonic interactions among phonons or from some other mechanisms 
such as the polariton coupling in ionic crystals!™ or polaron mechanism0. In such cases, 
equilibrium distribution of phonons are that of squeezed thermal phonons0. Therefore, it 
seems to be an important question how to determine the equilibrium distribution of phonons 
when there is a possibility that phonons can be found to be in non-classical states. As a 
particular example of some considerable interest, the squeezed states of phonons due to the 
photon - optical phonon interaction in an ionic crystalU should be mentioned here. The 
polariton coupling in such a system is described by the following Hamiltoniant3 

z k 

H k = uj k a[a k + u b b\b k + ig k [{a\ - a_ k ){b\ + 6_ fc ) + (at k - a k ){b_ k + &£)] 

where uj k is the photon frequency, u>b is the frequency of transversal oscillations of opti- 
cal phonons, g k is the polariton coupling constant and the operators a k ,b k describe the 
annihilation of photons and optical phonons respectively. Since the Hamiltonian under con- 
sideration is the Hermitian bilinear form, it can be diagonalized by the Bogolubov canonical 
transformation! similar to that used in the definition of squeezed statesa As a result, the 
thermal equilibrium state of the system is described by the following density matrix 



where H p denotes the Hamiltonian H in diagonal (polariton) representation and f3 is the 
reciprocal temperature. In analogy to the quantum optics, consider the so-called degree of 
coherenceEl 

G ( 2 ) = <W> 
' (6t 6 > 
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where (...) denotes the average with respect to the density matrix p(/3). It is straightforward 
to calculate as a function of temperature for typical parameters of an ionic crystal 
(see Fig. 1). One can see that, at low temperatures, ~ 8, while the same correlation 

(2) 

function calculated with the Bose-Einstein distribution gives 

G k b' e = 2- It is also seen that 

the strong quantum fluctuations can be observed only below T ~ 50A" because they are 
eroded by thermal fluctuations with the increase of temperature. 

In contrast to the case of non-classical states of photons there is no an efficient direct 
method of measurement allowing the characterization of the quantum state of Bose-type 
excitations in solidsB Even though correlation functions to any order would be demanded 
to describe fully a quantum state, it is usually good enough to distinguish quantum states 
by their number variances^!. Here, we present a way to determine the number variance of 
phonons at equilibrium in a Raman active medium. It is already suggested that correla- 
tion Raman spectroscopy may be used to measure the quantum statistical properties of a 
vibration mode for the case of Stokes (S) type Raman scatteringthrough a measurement 
of the intensity and the Mandel's Q-factor of the Rayleigh modeEl However, even at low 
temperatures vacuum fluctuations of the Anti-Stokes (AS) modes might disturb measure- 
ments of high order correlations and thus careful study of the role of the AS modes in such 
measurements is demanded. In this article, we follow a similar ideology in more general 
terms by examining both the S and AS components of multi-mode Raman scattering. Even 
though the problem becomes analytically intractable when AS modes are included, it is now 
possible to establish an interesting connection between the number variance of phonons and 
the correlations of S and AS modes. Moreover, due to the removing low temperature re- 
striction in the exclusion of AS modes, influence of temperature in the high order quantum 
correlations can be examined as well. 

The paper is outlined as follows. In Sec. II, using a general model of Raman type three- 
body scattering, we find the inter-mode correlation function of S and AS modes. Discussion 
of this general result under standard approximations of Raman scattering, with an emphasis 
of modifications in their range of validity, is subject to Sec. III. Finally, Sec. IV gives a brief 
summary of our results and conclusions. 

II. CORRELATION OF STOKES AND ANTI-STOKES PHOTONS 

General relations between the correlation function of S and AS modes and the number 
variance of phonons is developed in this section for the following Raman-type Hamiltonian, 

H = ^kA^A^A + (^kk'q a k'S a ki? a qV + ^kk'q a k'A a ki? a qV +H.C.), (1) 

kA kk'q 

where a^xi^ux) are the creation (annihilation) operators for the A-mode with momentum k 
and corresponding frequency cj kA . Here the mode index A = S, A, V, R stands for Stokes, 
Anti-Stokes, vibration and Rayleigh modes, respectively. As usually, the polarization labels 
are suppressed within the momentum symbols for the sake of notational simplicity. Coupling 
constants are denoted by M kk , q for the S-type scattering and M kk , q for the AS-type scat- 
tering. While writing this tri-linear bosonic Hamiltonian we assumed as usuallyEZl that the 
Raman scattering is observed under the condition lor,s,a ^ when the pair-wise creation 



3 



of radiation modes has quite small probability so that energy is conserved. This supposition 
is equivalent to the rotating wave approximation of the quantum optics0. We also assumed 
that the radiation consists of three R, S, and AS pulses which are well-separated on the 
frequency domain so that [a^x, aj^,] = ^kk'^AA'- If a single-mode strong coherent (classical) 
pumping is assumed, all one can expect is that the phase-matching conditions would have 
limited the number of active phonon modes to one. Nevertheless, it seems to be reasonable 
to consider the Raman scattering by an infinite Markoffian system of phonons0'0. In par- 
ticular, it permits oneself to take into account the broadening of S and AS* lines. The usual 
selection rules of Raman scattering, namely phase-matching or quasi-resonance conditions^, 
are not essential for the derivation of the general relations below. Therefore, the results given 
in this section are also valid in not so perfect Raman coupling situations which should be 
important in real materials. 

If we define the number operator n^x for the A-mode with momentum k as n^x = a^a^x, 
then the total number operator N x for A-mode becomes Nx = J2k n k\- Heisenberg equations 
of motion yield the conservation laws, also known as Manley-Rowe relations^, 

N S + N A + N R = C 1 , (2) 
N S -N A -N V = C 2 . 

Here constant operators Ci,C 2 are specified by the initial conditions. Similar relations can 
also be constructed for the scattering of photons of a monochromatic laser beam from a 
dispersionless optical phono iSi. Solving these equations for N$ and N A , the S and AS 
correlation function is found to be 

< N A ; N s >= \(V(C X ) - V(C 2 ) + V(N R ) - V(N V ) -2<C x -N R >-2< C 2 ; N v >), (3) 

where the correlation function < A; B > of two operators A, B is defined by 

< A;B >=< AB > - < A>< B > 

and hence variance of operator A is given by the self-correlation function V(A) =< A; A >. 
Here the averages < . > are with respect to the initial state since Heisenberg picture is used. 
It is natural to consider an initial state in which the S and AS modes are in their vacuum 
states, when we obtain, 

< N A (t);N s (t) > = ~(V(N R (0)) - V(N v (0)) + V(N R {t)) - V(N v (t)) 

- 2 < ^(0); N R (t) > -2 < AV(0); N v (t) >). (4) 

An operator A at time t is indicated by A(t) while initially by A(0). That equation connects 
the S and AS correlation function to the quantum statistical behavior of phonons and pump 
photons. 

Within conventional Raman theory Quantum properties of pump are usually neglected 
through the classical pump assumptionEJH This approximation introduces a time range 
to the problem during which changes in the pump intensity remains negligible. We can 
apply a similar approximation by assuming an intense laser pump with photons in coherent 
states and performing a mean field average over them in the above equations. Under this 
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assumption, the correlation function of the S and AS modes is related only to phonon 
statistics and the initial, known, number variance of the pump photons. However, time range 
of validity for the parametric approximation should be modified in our case. As we shall 
show in the subsequent section, statistical behavior of the pump might change significantly 
in shorter time than the occurance of a significant change in its intensity. Our purpose is to 
examine the equilibrium statistics of phonons determined by V(Nv(0)); therefore we need 
to express all time dependent terms on the right hand side of the Eq.4 in terms of initial 
operators to see any further relation between the S and AS correlation function and the 
equilibrium variance of phonons. For that aim we specify a model system and study its 
dynamics. 

We conclude this section by noting that a similar relation can be derived for the molecular 
Raman model, which is equivalent to the full bosonic Raman model under Holstein-Primakoff 
approximation in the case of low excitation density^. In that case, S and AS correlations 
depend on the quantum statistics of population distributions of the molecular energy levels. 



III. DISCUSSIONS FOR PARAMETRIC RAMAN MODEL 

In reality, coupling of one vibration mode to the pump beam for sufficiently long time 
of measurement is not an easy task. Therefore, in this section we investigate a Raman 
scattering in which coupling of pump photons to all phonon modes are allowed. We shall 
treat the pump as an intense coherent beam of photons and thus its state | ipn) in general 
is described by a multimode coherent state, 

I fa) = II ® I «i) ( 5 ) 
1 

in which ot\ are the coherence parameters of the modes 1. According to the remarks at the 
end of previous section, we now perform mean field averaging with respect to pump photon 
states in Eq.l assuming the Raman-active material is placed in an ideal cavity which selects 
single modes for S and AS radiations, namely k' = kA,s- Then after dropping constant terms 
the Hamiltonian in Eq.l reduces to an effective one, 

H eff _ J2 LU X n X + J2 ^qVaqV a qV + XXSq^qV + ^A^qY + H.C.), (6) 
X=S,A q q 

where new effective coupling constants g^' s are introduced by 

k 

The summation above can be calculated once the density of states for the pump is also 
specified. As one can see, the Hamiltonian will be in the given form, involving summations 
over phonon modes, in all cases except the case of perfectly phase matched single pump 
and phonon modes. In order to make sure that our results are not too susceptible to any 
imperfectness of the system arising from multi-mode nature of pump or phase-mismatches 
among the phonon and photon modes, we shall treat the problem using the model described 
by the above Hamiltonian involving summations over phonon modes. When finite number 
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of phonon modes are assumed, which is reasonable for real crystals of finite size, then such a 
model becomes integrable since the dynamics is ruled by the following closed set of operator 
linear differential equations, 



* j£ fl qv = ^qv« q v + 5^4 + 9q* a A, 
d 



! - ~UJ S al - X)#q* Q qV, (8) 



dt b q 

i ~d\ CLA = UaCLa + ^2 9q a civ- 

Let us introduce a vector of operators such that Y = [a' s , cla, { a qv}] T - We denote the matrix 
of coefficients in the above set of equations by M and its diagonalizing matrix by D, so that 
D~ l MD = El with eigenvalues E. Thus, we get 

Yi(t) = DijDj£Y h {0) exp (-lEjt), (9) 

where summation over repeated index is implied. It is therefore possible to write the solution 
for A = S, A-modes in the form, 

a A (t) f = u x (t)a j s + v x (t)a A + ^vOOaqV- ( 10 ) 

q 

Operators without time arguments are taken at t — 0. Time dependent parameters u, v, w 
are determined by the elements of matrix D and eigenvalues E. Let us note here that some 
general relations exists among u,v,w due to the commutation relations for a\ operators and 
they are not independent each other. More explicit way of evaluating u,v,w is presented 
below for the single mode phonon case where vector Y reduces to three dimensions in 
operator space. When there are no scattered light modes initially, the correlation function 
of S and AS modes becomes 

(n s {t);n A {t)) = A(t) + ^ 5 kq (t)(4 y a q y) + Cki Pq (t)(4v a qv; a\v a pv)- (H) 

kq klpq 

Here, parameters A,B,C are functions of u,v,w. Since the summations above can be 
converted into integrals involving phonon density of states, we see that if there are Van 
Hove singularities corresponding to the modes selected by Raman scattering, as in the case of 
recent experiments on the generation of non-classical phonon states via Raman scatterings^, 
then the correlation of S and AS modes will be determined strongly by that mode. If this 
not the case, then one can still expect domination of the modes obeying Raman selection 
rules. Then for that mode the random phase approximation permits us to write0 

(n s (t)} = \v s (t)\ 2 + \w' s \ 2 (l + n v ) (12) 
(n A (t)) = | u A (t) | 2 + | w' A | 2 n v 
(n s (t); n A {t)) = A'(t) + B\t)n v + C\t)V(n v ), 

in which the momentum label corresponding to relevant mode is fixed and dropped for 
the notational simplicity and primed parameters evaluated at that mode. It is possible 
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to argue by the results above that a measurement of the correlation between S and AS 
can be utilized to determine the variance of vibration modes, which we usually consider as 
phonons here, provided one knows the mean number of such modes initially. The latter 
information can be determined by either one of the first two relations in Eq. 12 after 
measurement of radiation mode intensities. Also measurement of radiation mode intensities 
and the knowledge of initial phonon number allow one to keep track of the evolution of mean 
phonon number through the Manley-Rowe relations given by the Eq.(2). Interestingly, since 
the mean number of phonons with non-classical distributions deviate significantly from that 
of Bose-Einstein distribution, it might be possible to find some traces of non-classicality 
even here. However, in order to classify the distribution of phonons strictly it would still 
be necessary to find the next moment of the distribution, in other words the variance of 
phonons. 

Now, an explicit way of determining u, v, w parameters will be demonstrated for the case 
of a single phonon mode. Because of three dimensional operator space in this situation, 
eigenvalues E t are found to be as the roots of the cubic equation 

£3 + 3uyE 2 _ {uJ 2 R _ 3u 2 y + ( |^|2 _ + + Uy) + \g A \\u R - U V )\ + 

+ UO V (uOy - uo 2 R ) = 0. 

Introducing coefficients PuQi as 

(Ei + uv)(Ei +UJR + uj v ) + |/| 2 - \g A \ 2 



Pi = 

Qi = 



g s Pi + E l + u v 



g A* 



we write the field operators as 

a^HEW**, (is) 
i 

a A {t) = Y.QiAie lElt . 
i 

Common operator coefficients A\ are determined in terms of the operators ay(0), a' s (0), o,a{0) 
using the Cramer's rule Ai = det (A)/ det (D) where 



D 



/111 

Pi P2 Ps 



and D\ is the matrix obtained by replacing the elements in the Zth column of D by the column 
vector [ay(0), a|(0), ^(O)] 7 . Thus, parameters u, v, w are determined in terms of interaction 
constants and the frequencies. More explicit expressions are too long and not very illumi- 
nating to reproduce here, but above analysis is quite suitable for numerical computation 
when some experimental data is available. At that moment we shall content ourselves with 
more fundamental discussions only. 
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In order to give a brief discussion of the dependence of the correlation function in Eq.9 
on squeezing parameter and temperature, we consider an equilibrium distribution of vibra- 
tion mode as of the squeezed thermal state with the following mean number and number 



variance^ 



( n v) — n v cosh 2r + sinh 2 r, 
Po(%) = {ny + ny) cosh4r + - sinh 2 2r, (14) 

where fiy is the mean number of phonons according to Bose-Einstein (BE) distribution 
and r is the real squeezing parameter. When r = 0, we recover the usual BE-distribution. 
According to Eq.10 the S and AS correlations increases with variance of phonons. And 
since both the riy and the V{ny) increases with temperature due to Eq.12, we see that 
temperature enforces stronger correlations of S and AS modes. However, we need to put a 
word of caution here, since the fluctuations which are determined by the self-correlations of 
the modes also increases with the temperature. In order to represent this competition, one 
can consider the cross-correlation function defined by@ 

t'S-AS = i ■ (15) 

y/V(n s )V(n A ) 

Since the denominator can be expressed in a similar structure as with the correlation func- 
tion in Eq.10, the cross correlation function will eventually saturate at high temperatures 
and at high squeezing parameters. Therefore, at high temperatures thermal fluctuations 
becomes important but not more important than in any typical quantum measurement. An 
estimation for a typical ionic crystal, for example, shows that the level of quantum fluctua- 
tions of phonon number exceeds that of thermal fluctuations below 30 -j- 50Kffi We also see 
through Eq.10 and Eq.12 that S and AS correlation increases with the squeezing parameter 
r. 

Finally, we examine the time range of validity for the parametric approximation. For 
that aim, we consider the Hamiltonian given in (1) for the case of perfect coupling of single 
modes. Let us suppress the momentum within the mode labels R, S, V, A and calculate a R (t) 
for times close to the beginning of interaction^. Up to the second order, we get 

a R (t) = e~ iu)Rt {a R + it(M s *a s a v + M A *a A a v ) - U 2 {\M s \ 2 v + \M A \ 2 ^)) (16) 

where v = a R (ns + ny + 1), \i = aR(n A — ny). Here, operators at t — are those without 
time arguments. Then, we calculate the mean number and the variance of pump photons 
for S and AS modes are in vacuum states initially as 

n R (t) =n R - t 2 (\M s \ 2 n R (l + n v ) + \M A \ 2 n R n s ), (17) 
V(n R (t)) = V(n R ) + 2t 2 {\M s \ 2 {V{n R ){\ + n v ) + n R (l + n v ) + 
+ \M A \ 2 {V{n R )n v - n R n v ))) 

In these equations averaging symbol, < . >, is not shown. Using the relation V(n R ) = n R 
for a coherent field, we find the time ranges t << ti,t 2 , for which the field intensity and the 
variance remain close to their initial values, as 



S 



1 \MS\i(l+n v ) + \M A \in v ' v ; 

1 

T2 = A\M s \ 2 (l + n v )' 

Clearly, we see a rescaling of time range of the usual time range of parametric approximation. 
At low temperatures ny ~ and thus r 2 = (1/4)^ shows a reduction of time range to 1/4 
of the typical range of parametric approximation. As an estimation, we may take g s ~ 
10 7 H&B, giving time ranges as T\ = 10 fs and r 2 = 2.5 fs. These ranges are readily available 
due to the remarkable recent developments in the field of femto-second spectroscopy©!^ 



IV. CONCLUSION 

Summing up our results, we should stress that the measurement of Stokes-anti-Stokes 
correlations looks like a reasonable method for detecting the number variance of a Raman- 
active vibration mode in solids. The most interesting and crucial fact is that the above 
method permits us to determine the number variance at thermal equilibrium, in other words, 
the variance just before the application of the pump beam. The phonon sub-system could 
be be in a non-classical state due to an interaction providing necessary correlations among 
phonons before the pump beam is applied. That interaction could be some anharmonic 
coupling with the heat bath, polaron or polariton mechanisms. Since these mechanisms 
are usually weaker than the first order Raman effect, after the application of the pump 
beam, dynamics of the phonon system is governed mainly by the Raman effect. Therefore, 
initial non-classical state of phonons and non-classical effects like squeezing which require 
phase coherence might be destroyed. That is why we have determined the general and 
fundamental formula given by Eq.4 in terms of the initial state of phonons and showed that 
under certain conditions it provides direct information on the initial, thermal equilibrium 
variance of phonons. Analyzing those conditions of applicability, we propose that at liquid N 
temperatures, using an intense coherent beam of ultra-fast laser source such as Ti-sapphire as 
a pump for a Raman active medium, one can measure the number correlation of the scattered 
Stokes and anti-Stokes modes and the mean photon numbers in these modes simultaneously 
by some photon counters, in order to determine the number variance of the vibration mode 
at equilibrium. The measurement can be realized through the use of a homodyne-type 
scheme^ in which the 5* and AS photons are counted by two different detectors connected 
with a computer fixing the simultaneous arrival of the 5* and AS photons. It is also shown 
that when the vibration mode is in squeezed state then an increase in the correlation of the 
Stokes and anti-Stokes modes occurs. 

Case of a multi-mode pump, important for ultra-short pulses, can be handled easily for 
materials which involves a strongly preferred phonon mode due to a Van Hove singularity 
in the frequency range of the pump, by an appropriate calculation of the effective coupling 
constants defined by Eq.7 which in turn modify only the coefficients A', B', C in Eq.12. Thus 
our conclusions should also be valid in this case. For materials in which such phonon modes 
are many or not exists at all, then application of a multi-mode pump and measurement 
of Stokes-Anti-Stokes correlation would still provide information on multi-mode phonon 
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correlations according to the general formula Eq.ll. This is a valuable knowledge to classify 
a possible non-classical multi-mode state of phonons like a multi-mode squeezed state. 

So far, the best achievement in squeezing of phonons is reported to be O.O1%0, provided 
by second order Raman scattering. We would like to emphasize that this is not the squeezing 
parameter r of Eq.12 but related to V(ny). Hence, the change in the Stokes-Anti-Stokes 
correlations we expect to be in the same order. There are other mechanisms which result 
in non-classical excitations in solids with different expressions and larger values for r and 
V(n). In fact, squeezing parameter reflects the strength of interaction preparing the non- 
classical state of these excitation! which is the initial phonon state in our scheme. The 
example of optical polariton we have discussed in the introduction, provides a two-mode 
squeezed state with squeezing parameter in the range r ~ 0.1 — 0.01 in Cud. Therefore, 
such a measurement with the ultrafast Raman correlation spectroscopy should not be too 
challenging and looks promising in our opinion. 

Let us finally note that the case of molecular Raman spectroscopy can also be treated 
with a similar formalism to get information on the quantum statistics of populations of 
molecular energy levels. 
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FIGURES 

FIG. 1. Phonon degree of coherence versus temperature for typical parameters of an ionic 
crystal: Q = 200K, g = 25K. 
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